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We consider the continuous measurement of a double quan- 
tum dot by a weakly coupled detector (tunnel point contact 
nearby). While the conventional approach describes the grad- 
ual system decoherence due to the measurement, we study the 
situation when the detector output is explicitly recorded that 
leads to the opposite effect: gradual purification of the double- 
dot density matrix. Nonlinear Langevin equation is derived 
for the random evolution of the density matrix which is re- 
flected and caused by the stochastic detector output. Grad- 
ual collapse, gradual puriflcation, and quantum Zeno effect 
are naturally described by the equation. We also discuss the 
possible experiments to confirm the theory. 



The problem of quantum measurements has a long 
history, however, it still attracts considerable attention 
and even causes some controversy, mainly concerning the 
wavefunction "collapse" (see, e.g.Jj'B). Among various 
modern approaches to this problem let us mention the 
idea of replacing the collapse postulate by the gradual 
decoherence of the density matrix due to the interaction 
with the detectom and the approach. oLa stochastic evo- 
lution of the wavefunction (see, e.g.,a~Eil). The latter ap- 
proach (which is used in the present paper) can describe 
the selective measurements for which the system evolu- 
tion is conditioned on the particular measurement result 
(other keywords of the approach are: quantum trajecto- 
ries, quantum state diffusion, quantum jumps, etc.). The 
renewed interest in the measurement problem is justified 
by the development of experimental technique, which al- 
lows more and more experimental studies of quantum 
measurement in optics and mesoscopic structures.E3~tII 
The problem also has a close connection to the rapidly 
growing fields of quantum cryptography and quantum 
computing.tS 

In the recent experimenttj with "which-path" interfer- 
ometer the suppression of Aharonov-Bohm interference 
due to the detection of which path an electron chooses, 
was observed. The weakly coupled quantum point con- 
tact was used as a detector. The interference suppression 
in this experiment can be quantitatively descriibsdJDy the 
decoherence due to the measurement processOcj 

We will consider a somewhat different setup: two quan- 
tum dots occupied by one electron and a weakly coupled 
detector (point contact nearby) measuring the position 
of the electron. The decoherence of the double-dot den- 
sity matrix due to continuous measurement in this setup 



has been analyzed in Refs.t2lL3. However, the decoher- 
ence approach cannot describe the detector output that 
is a separate problem analyzed in the present paper. We 
answer two interrelated questions: how the detector cur- 
rent behaves in time and what is the proper double-dot 
density matrix for a particular detector output. We pitcpw 
that the models of point contact considered in Refs.E^"Ell 
describe an ideal detector. In this case the density ma- 
trix decoherence is just a consequence of averaging over 
all possible measurement results. For any particular de- 
tector output our equations allow the evolution of pure 
wavefunction to be followed. Moreover, a mixed density 
matrix can be gradually purified in the course of a con- 
tinuous measurenierit. 

Similar to Ref.tj let us describe the double-dot sys- 
tem and the measuring point contact by the Hamiltonian 
T-L = Hdd + Ti-pc + T-Ltnt, where TLdd = (e/2)(c|ci - 
c\c2) + H(c\c2 + c\ci) is the Hamiltonian of the double- 
dot, Hpc = I]; Eia\ai+Y,r Eralo'r+Y.i,r T{alai+ajar) 
describes the tunneling through the point contact (T and 

H are real), and Hint = J2i r '^I'^^i'^t^i + o-jo'v), i-e. 
the tunneling matrix clement for the point contact is T 
or T + AT depending on which dot is occupied. So, 
the average current Ii = piprC^V /% flows through 
the detector when the electron is in the first dot {V is 
voltage across the tunnel contact, pi and pr are the den- 
sities of states), while the current is I2 — h + AI ~ 
2tt{T + AT)'^piPre'^V/h when the second dot is occupied. 

We make an important assumption of weak coupling 
between the double-dot and the detector (a better term 
would be the "weakly responding" detector), 

|A/| «/o = (/i+/2)/2, (1) 

so that many electrons, N > (/q/A/)^ ^ 1, should 
pass through the point contact before one can distin- 
guish which dot is occupied. This assumption allows the 
classical description of the detector, namely, to neglect 
the coherence between the quantum states with different 
number of electrons passed thro ugh t he detect or .cj 

The decoherence rate — (-^//i/e— ^/2/e)^/2 of the 
double-dot density matrix a(t) due to the measurement 
by tunnel point contact has been calculated in Ref.ta. 
In the weakly-responding limit (|l|) it can be replaced by 
Td — (A/)^/8e/o or by the expression 

- {Alf/ASi, (2) 

where Si — 2e/o is the usual Schottky formula for the 
detector shot noise spectral density Sj. Equation (0) has 
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also been obtained in RefsB3~H^ for the quantum point 
contact as a detector, the difference in that case is Si 
2e/o(l — T) where T is the transparency of tfie,channelcil 
(while above we implicitly assumed T <^ l).c§E3 Motice. 
that the decoherence rate (j|) was derived in Refs.t^TEHl 
without any account of the information provided by the 
detector, implicitly assuming that the measurement re- 
sult is just ignored. Now let us study how this additional 
information affects the double-dot density matrix. 

We start with the completely classical case in which 
there is no tunneling between dots {H = 0) and the initial 
density matrix of the system does not have nondiagonal 
elements, cri2(0) = (Ti2(i) = 0. We can assume that the 
electron is actually located in one of the dots, but it is not 
known in which one, and that is why we use probabilities 
CTii(O) and 0-22(0) = 1 — crii(O). The detector output is 
the fluctuating current I{t). The fluctuations grow when 
I{t) is examined at smaller time scales, so some averaging 
in time ( "low-pass filtering" ) is necessary, at least in or- 
der to neglect the problem of individual electrons passing 
through the point contact. Let us always work at suffi- 
ciently low frequencies, / ~ t~-^ <^ Si/e^, for which the 
possible frequency dependence of Sj can be neglected. 

The probability P to have a particular value for the 
current averaged over time r, (/) = JJ" I{t)dt, is 
given by the distribution 

P((/),r) =an(0)Pi((/),T) +(722(0) P2((/),r), (3) 
P.((/), r) = i27TD)-'/' exp [-((/) - I, f/2D] , (4) 

where D = Si/2t. Notice that these equations obvi- 
ously do not change if we divide the time interval r into 
pieces and integrate over all possible average currents for 
each piece (to consider only positive currents, the typ- 
ical timescale r should be sufficiently long, Sj/t <C P, 
that is always satisfied within the assumed low frequency 
range) . 

After the measurement during time r we acquire addi- 
tional knowledge about the system and should change 
the probj^bilities an according to the standard Bayes 
formulaEZI for a posteriori probability (taking into account 
particular detector result (/): 

aii(T)= (7n(0)exp[-((/) - /i)V2i?] 
x{aii(0)exph((/)-/i)V2i^] 

+ 022(0) exp[-((/)-/2)V2Z?]}'\ 

0-22(t)= f - O-ii(t). (5) 

Notice that we have considered so far the purely clas- 
sical measurement and did not use any "collapse" postu- 
late. Nevertheless, Eq. (^ can be interpreted as a grad- 
ual "localization" of the electron in one of the dots due 
to acquired information. 

Now let us ass ume that the initial state is fully coher- 
ent, |cti2(0)| = Vcrii(0)cr22(0), while stiU H ^ 0. Since 
the detector is sensitive only to the position of an elec- 
tron, the detector current will behave exactly the same 



way as in the case above. So, after the measurement dur- 
ing time r we should assign the same values for an (r) as 
in Eq. (||), but the question is not so trivial for the non- 
diagonal matrix element cti2(t). Nevertheless, we can 
write the upper bound: 
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(t)| < y/aii{T)a22iT). 



(6) 



If the actual measurement result is disregarded, then 
the upper bound for \ai2\ can be calculated using the 
probability distribution of different outcomes given by 
Eq. (^) and the upper bound (^ for each realization: 



I(^i2(r)). 



< 



v/^ll(T)a22(T)F((/),T)d(/) 



Vf^ii(0)fT22(0) exp[-(A/)V/45/] 



(7) 



(here ( )r means averagina-Qver realizations). This up- 
per bound exactly coincidesH with the result given by the 
decoherence approach (||). This fact forces us to accept 
the somewhat surprising statement that Eq. (^ gives not 
only the upper bound, but the true value of the nondi- 
agonal matrix element, i.e. the pure state remains pure 
(no decoherence occurs) during each particular measure- 
ment. (ActualW,|-tliis is the usual statement for selective 
measurementsjsilil i.e. when the detector output is taken 
into account.) 

Simultaneously, we proved that the^Qint contact de- 
tector considered theoretieally in Refs.EjEj (the model is 
confirmed experimentallylla) causes the slowest possible 
decoherence of the measured system, and hence repre- 
sents an ideaL detector in this sense. In contrast, tha 
result of RefH3 shows that a single-electron transistorEEl 
biased by relatively large voltage is not an ideal detec- 
tor Lihe non-ideal detector has also been considered in 
Ref.c3). Nptice, however, that in the range of elastic 
cotunneling^ the operation of the single-electron transis- 
tor is almost equivalentc^l to the case considered above, 
and, hence, it becomes an ideal detector. 

If the initial state of the double - dot is not purely co- 
herent, |fTi2(0)| < \/(Tii(0)cr22(0), it can be treated as 
the statistical combination of purely coherent and purely 
incoherent states with the same crii(O) and (722(0), then 



0-12 (r) = 0-12(0) exp 



oii(t)o22(t) 
011(0)022(0) 



1/2 



(8) 



Equations (||) and (||) are the central result of the 
present paper; they give the density matrix of the mea- 
sured system (in thfi-.case H = 0) with account of the 
measurement result .e3 

These equations can be also used to simulate the de- 
tector output I(t) and the corresponding evolution of the 
density matrix. For example, in the Monte-Carlo method 
we should first choose the timestep r satisfying inequal- 
ities e^/Si -C T ^ Si/(AI)'^ and draw a random num- 
ber for (/) according to the distribution (||). Then we 
update aii{t) and 022 (t) using this value of (/) and re- 
peat the procedure many times [the distribution for the 
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FIG. 1. Thick line: particular Monte-Carlo realiza- 
tion of (Jii evolution in time during the measurement of 
uncoupled dots, H = 0. Initial state is symmetric, 
cii(O) = 0-22(0) = 1/2, while the measurement leads to grad- 
ual localization. Initially pure wavefunction remains pure at 
any time t. Thin line shows the corresponding detector cur- 
rent (I) averaged over the whole time interval starting from 
t = while the dashed line is the current averaged over the 
running window with duration Si/{AI)^. 



current averaged over the interval At — t is new every 
timestep because of changing an (t) which are used in Eq. 
(^]. The nondiagonal matrix element can be calculated 
at any time with Eq. (||). 

Using Eqs. (||)-(||), this Monte-Carlo procedure can be 
easily reduced to the following nonlinear Langevin-type 
equation (equation for an is sufhcient): 



an =TZ, TZ 



2AI 



2AI 



-0-110-22 ■ 



Si 

<^22 - Cii 



Si 



Ai+m 



(9) 
(10) 



where the random process ^(t) has zero average and 
"white" spectral density S^ = Si. The second expres- 
sion for TZ allows the measurement to be simulated while 
the first one can be used to calculate the density matrix 
for given I{t) [in case _ff = it can more easily be done 
using Eq._p)]. 

Figure Hshows a particular result of the Monte-Carlo 
simulation for the symmetric initial state, o-ii(O) = 
o'22(0) = 1/2. The thick line shows the random evolution 
of CTii(i). Equation ( p^ ) describes the gradual localiza- 
tion in one of the dots (first dot in case of Fig. |^). Let us 
define the typical localization time as tioc — 25//(A/)^ 
(we choose the exponential factor at an = a^i = 1/2). 
Then it is exactly equal to the time Tdis = 25//(A/)^ 
necessary to distinguish between two states (defined as 
the shift of two Gaussians (^ from /q by one standard 
deviation), and tjoc = Td/2, where Td = FJ^. It is easy to 
prove that the probability of final localization in the first 
dot is equal to o-ii(O), because aaij) averaged over re- 
alizations is conserved (the deterministic flow of an due 
to the first term in square brackets of Eq. (HO) is exactly 



canceled on average by the dependence of the diffusion 
coefficient on an)- 

The detector current I{t) basically follows the evolu- 
tion of auit) but also contains the noise which depends 
on the bandwidth. The dashed line in Fig. [l] shows the 
current {I{t, t — At)) — At^^ It-At ^^^^ '^^ averaged over 
the "running window" with duration At — Si/{AI)'^, 
while the thin solid line is the current {I{t, 0)) averaged 
starting from t = 0. 

Now let us consider the general case of the double- 
dot with non-zero tunneling H. If the frequency f2 of 
"internal" oscillations is sufficiently low, = {4H^ + 
e^)^/^/7i <^ Si/e^, we can use the same formalism just 
adding the evolution due to finite H (the product ^tioc 
is arbitrary) . A particular realization can be either simu- 
lated by Monte-Carlo procedure similar to that outlined 
above [now update of ai2 {t) using Eq. (^) should be nec- 
essarily done at each timestep, together with the evo- 
lution due to finite H] or equivalently described by the 
coupled Langevin equations 



(Til = -0-22 = {-2H/h) Im((7i2) + TZ, 

ie iH , s 0-22 - an 



ill) 



h 

-7dO-12, 



— (CTii — 0-22) 



20-11022 



TZa 



12 



(12) 



where 7d = for an ideal detector (see below). The 
alternative "microscopic" derivation of these_£quations 
can be done for the particular model of Ref.Efl and will 
be presented elsewhere. 

Notice that in Eqs. (|9|)-([l^) the derivative is defined 
a.s a{t) — limT-^o[o'(i + T/2) — o(t — r/2)l/r (Stratonovich 
formulation of the stochastic equation^). The equations 
would be different if the definition [a{t + T) — a{t)]/T was 
used (Ito formulation). We use the former one because 
it gives the correct limit when the noise term ^{t) is re- 
placed by a sequence of smooth functiong23 and also be- 
cause the equations in Stratonovich formalism are physi- 
cally more transparent since they do not contain extra 
terms arising due to TZ^dt — const (for example, the 
usual calculus rule i^)' = /'<? + fg' is still valid). To 
translate Eqs. (||)~(|l2)-|into Ito formalism, one would 
need to add the termsEl {Si /2)F{dF/da)/2 where F is 
the factor before S,{t). This would lead to extra terms 
— (1722 — o'ii)A//2 in square brackets of Eqs. (p|)-(p^) 
and extra term -o-i2(A7)^/4S'/ in Eq. (|l|). Notice that 
in Ito formalism the equations become linear (except for 
the terms proportional to S,(t)). 

The simplest way to avoid the possible confusion be- 
tween two formulations of stochastic equations is to use 
the explicit calculation procedures (for finite t) described 
above. However, the difference should be taken into ac- 
count when results of otherj-approaches to the stochas- 
tic wavefunction evolutioru^OJ are compared. For exam- 
ple, this explains the apparent jxlifference between Eqs. 
(p|)-(p2[) and the results of RefHI for a two-level system 
(with £ = and 7^ = 0) derived in a different way. 
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FIG. 2. Random evolution of an (particular Monte-Carlo 
realizations) for asymmetric double-dot, e = H, with the elec- 
tron initially in the first dot, (Tii(0) = 1, for difi^erent strength 
of coupling with detector: C = h{AI)^ / Si H =0.3, 3, and 30 
from top to bottom. Dashed line represents C — (unmea- 
sured double-dot). Increasing coupling with detector destroys 
the quantum oscillations (while wavefunction remains pure at 
any t), slows down the transitions between states (Quantum 
Zeno effect), and for C ^ 1 leads to uncorrelated jumps be- 
tween well localized states. 



Among various approaches to selective quantum mea- 
surements, our approach is most clgsely related to the 
method of restricted path integraljlD however, in some 
sense we consider the classical (not quantum) path in- 
tegral. Let us also niention that the quantum nondemo- 
lition measurementaj are outside the scope of our study, 
we consider only the measurements at the so-called "stan- 
dard quantum limit" . 

Figure ^ shows the particular results of the Monte- 
Carlo simulations for the double-dot with e = H and 
the different strength of the interaction with an ideal de- 
tector. The electron is initially located in the first dot, 
o'ii(O) = 1- The dashed line shows the evolution of fxn 
without detector. Notice that because e 0, the ini- 
tial asymmetry of the electron location remains in this 
case for infinite time. When the interaction with de- 
tector, C — h{AT)'^/SiH, is relatively small (top sohd 
line), the evolution of cth is close to that without the 
detector. However, the electron gradually "forgets" the 
initial asymmetry and the evolution can be described as 
the slow variation of the two-parametric phase of oscilla- 
tions (recall that the wavefunction remains pure). In the 
decoherence approach (averaging ov£i; realizations) this 
corresponds to cth ^ 1/2 at i — > oo.Ej 

When the couphng with the detector increases, the 
evolution significantly changes (middle and bottom 
curves in Fig. ||). First, the tpjO^^timiipl^ dots slows 
down (Quantum Zeno effectH]a#ffl'BE3t3'El) . Second, 
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FIG. 3. Gradual purification of the density matrix [aii(t) 
and Imcri2 (t) are shown] of the symmetric double-dot (e — 0) 
measured by slightly nonideal {'jd/^d = 0.1) weakly coupled 
(C = 0.1) detector. 



while the frequency of transitions decreases with increas- 
ing interaction with the detector, the time of a transition 
also decreases, so eventually we can talk about uncorre- 
lated "quantum jumps" between states. 

In a regime of small coupling with a detector, C ^ 1, 
the detector output is too noisy to follow the evolution of 
an and, correspondingly, only slightly affects the oscilla- 
tions (the presence of quantum oscillations in the double- 
dot can be noticed only as a relatively small peak in the 
spectral density of the detector current). In contrast, 
when C ^ 1 the detector accurately indicates the posi- 
tion of electron and simultaneously destroys the oscilla- 
tions. 

Equations (pT)-(p^ with the term TZ given by Eq. (||) 
can be used to obtain the evolution of the density matrix 
in an experiment provided the detector output I(t) and 
initial condition (0) are known. Notice that even if 
the initial state is completely random, an = aii = 1/2, 
ai2 — 0, the nondiagonal matrix element gradually ap- 
pears during the measurement, so that sufficiently long 
observation with an ideal detector leads to almost pure 
wavefunction for the double-dot. Such a purification of 
the density matrix described by Eqs. (|l^)-(|l^) is analo- 
gous to the localization a,t H — 0. 

Equations (pT|)-(p^ can be generalizedjfor a nonideal 
detector, F^ >{Aiy/4Si (as in Refs£3'Ej), which gives 
less information than possible in principle. Let us model 
it as two ideal detectors "in parallel" with unaccessible 
output of the second detector. Then the information loss 
can be represented by the extra decohererice term — jdCu 
in Eq. (|l|) where 7^ = F^ - {AI)^/ASilB The hmiting 
case of a nonideal detector is the detector with no out- 
put (just an environment. A/ = 0) or with disregarded 
output. Then Eqs. (pT|)-(p^ reduce to the standard de- 
coherence approach. 

For a nonideal detector it is meaningful to keep our 
old definition of the localization time, tiqc — T^is — 
2Si /{AI)^ , while < 2tioc- So, we consider localiza- 
tion time not as a real physical quantity but as a quan- 
tity related to the observer's information. Similarly, the 
effective decoherence time is defined as = 7^^. 
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Figure ^ shows a particular realization of random evo- 
lution of CTii and Imai2 for a symmetric double dot mea- 
sured by weakly coupled (C — 0.1) nonideal detector with 
Id/^d = 0.1. We start from maximally mixed state, 
CTii(O) = 0.5, cri2(0) = 0, and the Figure shows the grad- 
ual purification of the density matrix in a course of mea- 
surement (notice that Re(Ti2(t) = because £ = 0). The 
nonideality of the detector does not allow the complete 
purification: oscillations of lm.a\2{t) do not reach ±0.5 
limit, as it would be in the case of ideal detector. 

Let us mention that following the "orthodox" (Copen- 
hagen) point of view, we do not attempt to distinguish 
between "real" density matrix and the density matrix 
which can be known by the observer. For example, the 
evolution of an due to the measurement in case of no 
tunneling between dots {H = 0) can be interpreted both 
as a real process or just as gradual acquiring of infor- 
mation about the electron position. Another example 
is the case of nonideal detector. We can interpret the 
term — 7^0-12 in Eq. (^ as real decoherence, however, it 
is also possible to argue that it just represents the par- 
tial loss of information inside the imperfect detector, so 
that perhaps the pure density matrix could be restored 
if some hidden traces left in the detector had been an- 
alyzed. Developing this example further, let us imagine 
that two observers have different levels of access to the 
detector information, then the density matrix for them 
will be different. Actually, this just means that the ob- 
server with less information will not be able to make as 
many (or as accurate) predictions as the other one. Nev- 
ertheless, he still can treat his density matrix as a real 
one for all purposes. The limiting case when the ob- 
server does not have any information about the detector 
output (or this information is ignored in the experiment) 
is equivalent to averaging over all possible realizations, 
i.e. to the standard decoherence approach. 

So, if different realizations of the detector outapt are 
effectively averaged in an experiment (as in RefJl3), the 
decoherence approach is suitable. In contrast, if the sin- 
gle realization of the detector current is recorded (and 
somehow used) in an experiment, then the proper de- 
scription is given by Eqs. (ll)-(p^). The simplest exper- 
imental idea is just to measure I{t) when C is not too 
small and check if it is consistent with these equations. 
However, it would be much more interesting to devise 
an experiment in which the subsequent system evolution 
depends on the preceding measurement result. 

For example, let us first prepare the double-dot in the 
symmetric coherent state, an — a22 — |o'i2| — 1/2, 
make H = (raise the barrier), and begin measure- 
ment with an almost ideal detector. According to our 
formalism, after some time r (the most interesting case 
is T ~ T/oc) the wavefunction remains pure but becomes 
asymmetric and can be calculated with Eqs. ^ and (H). 
To prove this, an experimentalist can use the knowl- 
edge of the wavefunction to move the electron into the 
first dot with probability equal to unity. Namely, he 
switches off the detector at i = r, reduces the bar- 



rier (to create finite H), and creates the energy differ- 
ence e = [(1 - 4|cri2p)i/2 _ l]i/Recri2/|cri2p; then after 
the time period At = [tt — arcsin(Imcri2 ?iri/-ff)]/n the 
"whole" electron will be moved to the first dot, that can 
be checked by the detector switched on again. (Alterna- 
tively, using the knowledge of aij (r) an experimentalist 
can produce exactly the ground state of the double-dot 
system and check it, for example, by photon absorption.) 

Another experimental idea is to demonstrate the grad- 
ual purification of the double-dot density matrix. Let us 
start with a completely random state {an = (T22 = 1/2, 
(^12 = 0) of the double-dot with finite H. Then using 
the detector output I{t) and Eqs. (|ll|)-(|l2|) it is possible 
to calculate the evolution of the density matrix. These 
calculations will show the gradual purification (the most 
interesting case is flrioc ^1), eventually ending up with 
almost pure wavefunction with precisely known phase of 
quantum oscillations. The final check of the wavefunc- 
tion can be similar to that considered above. However, it 
can be even simpler, because with the knowledge of the 
phase of oscillations it is easy to stop the evolution by 
raising the barrier when the electron is with certainty in 
the first dot. If rapid calculations (by some analog on- 
chip circuit) are not available, the barrier control can be 
random, while appropriate cases can be selected later. 

An experiment of this kind could verify the formalism 
developed in the present paper. While such an experi- 
ment is still a challenge for present-day technology, we 
hope that it can be realized in the near future. 

In conclusion, we have developed a simple formalism 
for the evolution of the double-dot density matrix with 
an account of the result of the continuous measurement 
by weakly coupled (weakly responding) point contact. In 
contrast to most previous studies on the selective quan- 
tum measurements, our equations treat mixed states and 
allow the consideration of a nonideal detector. The equa- 
tions show the gradual purification of initially mixed 
state of the double-dot due to continuous quantum mea- 
surement. This effect can be studied experimentally in 
various mesoscopic setups. 

The author thanks S. A. Gurvitz, K. K. Likharev, D. 
V. Averin, and T. V. Filippov for fruitful discussions. 
The work was supported in part by French MENRT 
(PAST), US AFOSR, and Russian RFBR. 
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